Regimes of classical transport of cold gases in a two-dimensional anisotropic disorder 
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We numerically study the dynamics of cold atoms in a two-dimensional disordered potential. We 
consider an anisotropic speckle potential and focus on the classical dynamics, which is relevant 
to some recent experiments. Firstly, we study the behavior of particles with a fixed energy and 
identify different transport regimes. For low energy, the particles are classically localized due to the 
absence of a percolating cluster. For high energy, the particles undergo normal diffusion and we 
show that the diffusion coefficients scale algebraically with the particle energy, with an anisotropy 
factor which significantly differs from that of the disordered potential. For intermediate energy, 
we find a transient sub-diffusive regime, which is relevant to the time scale of typical experiments. 
Secondly, we study the behavior of a cold-atomic gas with an arbitrary energy distribution, using 
the above results as a groundwork. We show that the density profile of the atomic cloud in the 
diffusion regime is strongly peaked and, in particular, that it is not Gaussian. Its behavior at large 
distances allows us to extract the energy-dependent diffusion coefficients from experimental density 
distributions. For a thermal cloud released into the disordered potential, we show that our numerical 
predictions are in agreement with experimental findings. Not only does this work give insights to 
recent experimental results, but it may also serve interpretation of future experiments searching for 
deviation from classical diffusion and traces of Anderson localization. 
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I. INTRODUCTION 

Transport processes are widespread and common in 
many fields of physics, chemistry and biology and have 
drawn the attention of both theorists and experimental- 
ists. Diffusion of waves and particles in disordered media 
determines the behaviors of a variety of physical systems. 
These include electronic conductivity in dirty materials 
PQ, superfluidity of liquid helium in porous substrates 
[2], and diffusion of light in dense media [3J, with ap- 
plications to dielectric materials [I] or interstellar clouds 
[5]. Diffusion is at the heart of transport phenomena 
in most systems of condensed-matter physics. For in- 
stance, it underlies the Drude theory of conductivity [B], 
as well as the self-consistent theory of Anderson localiza- 
tion [7J. A major hindrance to the complete understand- 
ing of condensed-matter systems is the complex interplay 
of many relevant ingredients, such as the structure and 
thermal fluctuations of the underlying substrate [SJ, the 
disorder [8], and the inter-particle interactions, which can 
lead either to superconductivity [S] or to metal-insulator 
transitions jTU], [TTJ . This poses a number of fundamen- 
tal questions regarding the existence and the nature of 
disorder-induced metal-insulator transitions |12j , the role 
of the spatial dimension and the effect of non-linearities. 

Currently, much attention is devoted to disordered 
quantum gases [T31ll7j . These systems offer unprece- 
dented control of most parameters, including tempera- 
ture, geometry and particle-particle interaction [T8ll21| . 
Moreover, disordered potentials can be engineered al- 
most at will, and their statistical properties can be con- 
trolled and tuned [22 . So far, a large body of the- 
oretical work has been devoted to Anderson localiza- 
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tion of non- interacting particles [23-35 , interacting Bose 
[36H47] and Fermi [48H5T] disordered gases, and two- 
component disordered systems [5"2"H5"8] . Experimentally, 
classical suppression of transport of Bose-Einstein con- 
densates [SHI E0] and Anderson localization of matter- 
waves [611 [62] in one dimension (ID) have been reported 
recently. The competition between interaction and dis- 
order has been investigated in ID bi-chromatic optical 
lattices [S31EI] and in three dimensional (3D) disordered 
lattices [651166]. 

Studying the dynamics of even non-interacting atoms 
in a disordered potential is a difficult problem, especially 
in dimensions higher than one. For classical particles, 
the motion in disordered potentials is generally chaotic, 
which sparks a variety of dynamical regimes [67] , includ- 
ing classical localization, normal diffusion as well as sub- 
and super-diffusion. For waves, diffusion is obtained in 
the limit of incoherent multiple scatterings 7 O.v 7oll7(ij . 
Multiple coherent scatterings [58H70] from the random 
defects can induce weak [72] or strong [5U] localization 
effects. With a view towards search for localization ef- 
fects in cold-atomic gases, the characterization of classi- 
cal transport regimes is a central task. It is of special 
interest in dimension two (2D), which is the marginal di- 
mension for return probability in Brownian motion and 
for Anderson localization [71]. In the context of cold- 
atomic gases, diffusion has been experimentally studied 
in random dissipative 3D optical molasses [731 El] j an d 
conservative 2D random potentials [77]. 

In this work, we numerically study the transport prop- 
erties of cold atoms in 2D disorder. We consider a speckle 
potential and focus on the classical regime such that 

A dB < o-r < Ib < L < L loc , (I) 

where er R is the characteristic length scale of the disor- 
dered potential, A dB is the atomic de Broglie wavelength, 
l B is the Boltzmann (transport) mean free path, L is 
the system size and L\ oc is the localization length. The 
inequality A dB <C l B guarantees that wave interference ef- 
fects, which are important in 2D, occur on a very large 
scale which will be assumed here to be much larger than 
the system size and thus disregarded. In particular, be- 
sides that, in 2D, all quantum states are localized [7T] . 
the inequality L <C L\ oc allows us to neglect strong lo- 
calization effects. When neglecting wave interference ef- 
fects, the main dynamical regime is provided by diffusion 
[3[7]. For er R < Z B (which is most often true), the con- 
dition a B < A dB may be fulfilled. In this case, the wave 
nature of the particles governs the scattering from each 
defect. Here instead we consider the opposite condition 
A dB <C cr R meaning that the atoms scatter in the disor- 
dered potential as purely classical particles, hence satis- 
fying the Newton equations of motion. Finally, since l B is 
the typical length of ballistic trajectories, the inequality 
l B -C L ensures that the number of scattering events in 
the system is sufficiently large to strongly affect the mo- 
tion of particles. As discussed in our paper, the regime 
indicated by Eq. (nj) is relevant to some recent cold-atom 



experiments [77] , 

The central aim of our paper is to analyze the trans- 
port regimes of classical particles at different energy. In 
particular, we show that they strongly depend on the 
topographic properties of the disordered potential. An 
additional aim is to use these results to determine the 
behavior of atomic clouds with a broad energy distribu- 
tion, as relevant for cold-atom experiments. 



A. Main results presented in the paper. 

Here, we give a general overview of the main results 
of this manuscript (leaving details for the text). Fig- 
ure [l] illustrates the results obtained for an anisotropic 
speckle potential of anisotropy factor A = a x /a y , with 
a r the correlation length in the r € {x, y} direction. At 
any energy, a particle first undergoes a ballistic flight for 
a typical time r* ~ Ib/v where v is its initial velocity. 
At sufficiently small energy, E, all classically allowed re- 
gions in the disordered potential have a finite size [see 
Fig. [TJa)], and the particle remains classically localized 
in the one containing its initial position. For a particle 
initially at the origin, the disorder-averaged mean square 
displacement, (r(t) 2 ) cx t ar , shows a sub-diffusive behav- 
ior (a r < 1) at intermediate times, which corresponds 
to the random exploration by the particle of its classi- 
cally allowed region. Once the full region is explored, 
(r(t) 2 ) saturates asymptotically in time (a r = 0) to a 
finite value [see Fig.[ljd)]. Being of purely topographical 
nature, this behavior leads to (x 2 )/(y 2 ) = A 2 . A critical 
energy E c signals the percolation phase transition corre- 
sponding to the appearance of infinitely extended allowed 
regions. For E slightly above E c , the percolating cluster 
(classically-allowed region) is characterized by large lakes 
connected by narrow bottlenecks [see Fig. gb)], which 
are located at the saddle points of the disordered po- 
tential. The particle then spends long times exploring a 
lake, with a similar dynamics as described above, lead- 
ing to the sub-diffusive behavior observed in Fig. [TJe) 
at intermediate (but experimentally relevant) times. For 
longer times and larger distances, the dynamics is domi- 
nated by rare but possible transitions between the lakes, 
which leads, asymptotically in time, to normal diffusion, 
(r(t) 2 ) — 2D r t. At larger energy, the lakes merge and the 
size of forbidden regions rapidly decreases [see Fig. life)]. 
Then, the dynamics is characterized by normal diffusion 
for t > t* [see Fig. [TJf )] . In this regime, we identify a 
power-law behavior of the diffusion coefficients as a func- 
tion of the particle energy and show that the anisotropy 
of diffusion significantly differs from the anisotropy of the 
disorder, D x /D y =^ A 2 . 

It is worth noting that a sub-diffusion regime should 
not be viewed as a precursor of localization. In fact, we 
see in the example considered here that sub-diffusion can 
cross over to either localization (for E < E c ) or normal 
diffusion (for E > E c ). 
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Figure 1: (color online) Topography and transport regimes of classical particles in a 2D, anisotropic, blue-detuned speckle 
potential with the anisotropy factor A = a x /a y = 2 (see Sec. II A I . Panels (a)-(c) show the allowed (white; V(r) < E) and 
forbidden (black; ^(r) > E) regions for a classical particle of energy E in the disordered potential V(r). Different columns 
correspond to the same realization of the disordered potential and different energies: (a) and (d) E = 0.3Va; (b) and (e) 
E = IVr; (c) and (f) E = 3Vr, where Vr is the average amplitude of the disordered potential (see Sec. II A I . In pane l (a) , 
the energy is below the percolation threshold E c — 0.52Vr (for the 2D blue-detuned speckle potential used here, see Sec. II B I 
and all the particles are classically localized. In panels (b) and (c), the energy is above the percolation threshold and there 
appear allowed regions connecting the opposite sides of the box. Panels (d)-(f) show the mean-square displacement, (r{t) 2 } 
(the direction r = x is represented by the upper thick blue line, r = y by th e lower thick red line) as a function of time in 
units of r R = ^/ma^/lVn], where m is the mass of the atoms (see Sec. Ill Al. The bl ack dotted line is the isotropic ballistic 
behavior, (r(t) 2 } oc t 2 predicted by Eq. ( 14 1. The dashed lines are (r(t) 2 } oc t r [see Eq. ( 13 1] for intermediate time scales giving 
a sub-diffusive behavior (a r < 1) in (d) and (e). The solid lines are {r(t) 2 ) oc t c " r for long times giving a time independent 
behavior (a r = 0) in (d) and a diffusive scaling (a r = 1) in (e) and (f). 



B. Outlook. 

The manuscript is organized as follows. In Sec. [TTJ 
we acquaint with the main statistical and topographic 
properties of the 2D speckle potential used along the pa- 
per. We identify a percolation phase transition and deter- 
mine the critical exponent characterizing the divergence 
of the classical localization length. In Sec. |III| we an- 
alyze the dynamics of non-interacting classical particles 
of given energy in the disordered potential. Firstly, we 
write down the classical equations of motion and, with 
a proper rescaling, we highlight the universality class of 
the classical dynamics. Due to the redistribution of ki- 
netic and potential energy induced by the disordered po- 
tential, the total particle energy is identified as the sole 
relevant parameter for the dynamics. Different transport 
regimes are then identified and discussed: (i) classical lo- 
calization, (ii) normal diffusion and (iii) transient subdif- 
fusion. Secondly, in Sec. |IV[ we use the above results as a 
groundwork to study the overall dynamics of atomic gases 
with a broad energy distribution, as relevant to cold- 



atom experiments. In particular, we show that, in the 
diffusion regime, the mean square size of the gas grows 
linearly in time ((r(t) 2 ) oc t) but, due to the summation 
of all diffusive components with energy-dependent diffu- 
sion coefficients, expanding density profiles significantly 
deviate from Gaussian functions. In Sec. |VJ considering 
a thermal cloud released in the disordered potential, we 
compare our predictions with the experimental results 
of Ref. [77]. We provide fitting functions to extract the 
single-energy diffusion coefficients from the collective ex- 
pansion of an atomic cloud of atoms. We find fair agree- 
ment between numerical calculations and experimental 
measurements. Finally, we summarize our findings and 
discuss possible extensions of this work in Sec. |VI| On 
the one hand, our analysis provides a guideline for the 
experimental investigation of classical transport in cold 
atomic gas in disordered potentials [77] . On the other 
hand, it may contribute to the interpretation of on-going 
experiments searching for deviation from classical diffu- 
sion and traces of Anderson localization with ultracold 
gases in 2D geometry. 
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Figure 2: (color online) (a) Example of anisotropic speckle 
potential used in the numerical simulations. The anisotropy 
factor is A = a x /a y — 2. (b) Autocorrelation function of 
the disordered potential (dots), along the x and y directions 
[C(x, y = 0) and C(x — 0,y), respectively]. It has been 
obtained numerically by averaging over 100000 realizations 
of the disordered potential. The solid lines correspond to 
Eq. (J3j) with no fitting parameters, (c) Intensity distribution 
obtained for a single realization of the disordered speckle po- 
tential (dots). The solid line is Eq. {2}. 



II. DISORDERED POTENTIAL 
A. Statistical properties of the speckle potential 

Throughout the paper, we consider a speckle disorder, 
as commonly devised in cold-atom experiments |22l 1591 - 
IBTHTTHSP] . A speckle pattern is created, for instance, by 
a coherent monochromatic laser beam passing through a 
diffusive plate [55] . A fully developed speckle field can be 
produced by a ground glass whose surface contains ran- 
dom grains associated to optical path length fluctuations 
longer than the laser wavelength. In this case, the speckle 
field, observed at the focal plane of a converging lens, re- 
sults from the coherent superposition of waves originating 
from different scattering sites on the rough surface with 
uncorrelated phases uniformly distributed in [0, 2ir]. Ow- 
ing to the central limit theorem, the real and imaginary 
parts of the electric field are independent Gaussian vari- 
ables. The light-shift potential V(r) felt by the atoms 
exposed to the electromagnetic radiation is proportional 
to the field intensity: it is thus a random potential but 
we stress that its statistics is not Gaussian. The speckle 
potential is also inversely proportional to the detuning of 
the laser light with respect to the two-level atomic transi- 
tion and can be either attractive or repulsive (see below) . 
Finally, assuming a sufficiently far-detuned laser field, we 
will neglect, in the following, dissipative effects. 

The statistical properties of the speckle potential have 
been extensively investigated [5T]. The probability that 
the potential, at a given point r, has amplitude V = V(r) 
follows an exponential law 



P{V) 



e~ v / v * (V 
\V R \ \V R 



(2) 



(V(r)) - V R and (AV(r)) 2 = (V(r) 2 ) - (V(r)) 2 = V 2 , 
where (...) indicates the average over the different realiza- 
tions of the disordered potential. The quantity V R is the 
amplitude of the disorder, and its sign depends on the de- 
tuning of the laser light from the atomic resonance: For a 
blue-detuned speckle potential, the disorder is repulsive 
(i.e. positive V R > and V(r) > 0) and it is characterized 
by arbitrary large intensity peaks with exponentially low 
probability. Conversely, for a red-detuned speckle poten- 
tial, the disorder is attractive (i.e. negative V R < and 
^( r ) < 0) and it is characterized by arbitrary low dips 
with exponentially low probability. 

The autocorrelation function of the disordered poten- 
tial, C(r) = (V(r+r')V(r')) - {V{r')) 2 , can be controlled 
by the geometry of the diffusive plate [S3] . Without loss 
of generality, it can be written as C(r) = V R c(r/a R ) |82j . 
where a R is the correlation length of the disordered po- 
tential, that is the typical width of C(r), which is pro- 
portional to the speckle grain size. The precise definition 
of <j r depends on the specific model of disorder. In this 
manuscript, we will focus on a 2D, anisotropic, speckle 
potential with a Gaussian correlation function of widths 
a x and a y along the two main directions, 



C(r) = Vl exp ( - x 2 /2a 2 x ) exp ( - y 2 /2a 



(3) 



We define a R = a x and A = a x /a y the anisotropy factor. 
In Fig. [2ja), we show a typical realization of a speckle 
potential as used in the numerical simulations (A = 2). 
Figures [2jb) and [2|c) show its autocorrelation function 
and intensity distribution, respectively. Experimentally, 
a speckle potential with a Gaussian correlation function 
can be obtained by illuminating a diffusive plate with a 
Gaussian laser beam. The anisotropy can be controlled 
by the aperture function of the diffusive plate [22] |8T] or 
by tilting a 2D, isotropic speckle field with respect to the 
diffusion plane of the atoms [77] . 



B. Topographic properties of the speckle potential 

We now study some topographic properties of the 2D 
speckle potential which are relevant to understand the 

Some 



transport regimes of classical particles [see Sec. Ill 



where 8 is the Heaviside step function. We have 



of them, as, for instance, the percolation threshold, which 
will be discussed below, do not depend on the anisotropy 
of the potential. Therefore, these are straightforwardly 
obtained from those of the isotropic case, upon rescaling 
of the y direction by the anisotropy factor A. We thus 
focus, without loss of generality, on an isotropic speckle 
disorder in this sub-section. 

The positivity of the kinetic energy constraints the tra- 
jectories of classical particles with total energy E to the 
sub-space defined by V(r) < E. This condition separates 
the space in allowed and forbidden regions. For the sake 
of illustration, Figs. [TJa)-(c) show the classically-allowed 
(V(r) < E; white) and forbidden (V(r) > E; black) 
regions for particles of fixed energy in a specific realiza- 
tion of a 2D blue-detuned anisotropic speckle potential. 
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In a continuous disordered potential, the topography of 
classically allowed regions strongly depends on the ratio 
2?/|V^|, and, for a speckle potential, on the sign of V R . 
The fraction of space satisfying the condition V(r) < E 
is given by [53] 

<t>{E) - [ dVP(V). (4) 

J — oo 

For a speckle potential, P(V) is given by Eq. (pi). For 
blue-detuned speckle potentials, we thus find 4>{E) = 1 — 
e -E/\v R \ for E > o and = (i.e. there is no allowed 
region) for _E < 0. For red-detuned speckle potentials, 
we have <f>(E) = e +E /\ VR \ for E < and = 1 (i.e. 

there are no forbidden region) for E > 0. 

For 2D disordered potentials, there exists a critical en- 
ergy, the so-called percolation threshold, E c , which marks 
the sharp transition between the existence of infinitely 
extended allowed and forbidden regions [83] . The perco- 
lation transition is often considered as the simplest possi- 
ble phase transition with nontrivial critical behavior |84j . 
It has a pure geometrical nature and exhibits long-range 
correlation near the critical value, corresponding to the 
divergence of the average size of the allowed regions (in 
infinite space). 

Below E c , the allowed regions are disconnected and 
form "lakes" of finite size [see Fig. Wa)]. The allowed 
regions consist of a finite number of minima of the po- 
tential which are connected through saddle points. In 
other words, particles of energy E < E c are classically 
localized, i.e. they cannot move to arbitrary large dis- 
tances. On a macroscopic scale, the system behaves as 
an insulator. Above E c , an infinite number of "lakes" 
merge together to form an infinite "ocean" : The allowed 
subspace contains at least one connected subset (cluster) 
which stretches to infinity [see Figs.[lpt>)-(c)]. In other 
words, particles with energy E > E c and initially placed 
in an infinitely-extended allowed region can move to in- 
finity. The system thus behaves as a conductor. 

A challenging issue is to determine the value of the 
percolation threshold E c for different models of disorder. 
In the case of a disordered potential with sign symme- 
try, [i.e. when the statistical properties of the potential 
V(r) are equivalent to those of 2(V(r)) — V(r)], the per- 
colation threshold is E c — (V(r)) and the corresponding 
critical area fraction is <f> c = 4>(E C ) — 1/2 [83 . An exam- 
ple is given by Gaussian random distributions. However, 
a speckle potential is non-Gaussian and does not have 
the sign symmetry discussed above. For instance, a blue- 
detuned speckle potential is bounded below by zero, but 
it is unbounded above (and vice versa for the red-detuned 
speckle potential) . We have numerically investigated the 
percolation threshold in 2D speckle potentials with the 
correlation function given by Eq. Q with A = 1. In 
particular, we have evaluated the probability (for differ- 
ent realizations of the disorder), P pC rc, to have an al- 
lowed region connecting the opposite sides of a square 
box. The results are shown in Fig. [3] for a blue-detuned 
speckle and for different box lengths, L. As expected, 
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Figure 3: (color online) Percolation probability as a function 
of the particle energy. It is calculated numerically as the 
probability to find at least one allowed region connecting the 
opposite sides of a square box of size L in a blue-detuned 
(isotropic) speckle potential. Different symbols correspond to 
different values of L (indicated in the figure) and the solid lines 
are guides to the eye. The vertical dashed line indicates the 
percolation threshold, E c — 0.52V R . The percolation proba- 
bility is calculated over 100 realizations of the disorder and 
the numerical grid step is o\r/10. 

the crossover between P per c = and P per c = 1 becomes 
sharper and sharper by increasing L, and all curves cross 
the value P pcr c = 0.5 at approximately the same value 
of E/V R . This is compatible with a phase transition 
at E c = 0.52Vr, for a blue-detuned, 2D, speckle poten- 
tial, in the limit L — > oo (vertical dashed line in the 
figure). According to Eq. Q, the critical fraction of 
allowed regions is given by <j) c = 0.405, a value signifi- 
cantly different from 0.5 as expected for potentials with 
sign-symmetry. The values E c = 0.52V R and 4> c = 0.405 
found here are in agreement with the experimental find- 
ings [55] and numerical results [55] obtained by mapping 
the speckle potential onto a network connecting potential 
minima through saddle points. Because of the duality 
between allowed and forbidden regions in the 2D case, 
when inverting the sign of the potential, allowed regions 
for a particle of energy E becomes forbidden regions for 
a particle of energy — E. The percolation threshold for 
the red-detuned speckle potential can then be directly 
calculated from that of the blue-detuned speckle poten- 
tial: E c {— V R ) — —E C (V R ). For red-detuned speckle po- 
tentials, we thus find that the percolation threshold is at 
E c = — 0.52\V R \ and the corresponding fraction of allowed 
regions is C = 0.595. 

In order to further characterize the topographic prop- 
erties of the speckle potential in a homogeneous box we 
have studied the classical localization length Jz?(E), i.e. 
the mean diameter of the percolating cluster stretching 
from the box center. We use the definition 

Se{E) = 2V^2, (5) 

where 5r 2 = J d 2 r A E (r)\r\ 2 / J d 2 r A E (r) with A E (r) 
1 in the classically allowed region stretching from the 
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Figure 4: (color online) (a) Probability distribution of 
the classical localization length denned in Eq. (Tsh , for a 
blue-detuned speckle potential. Symbols are numerical re- 
sults for different energies (indicated in the figure). The 
solid lines are fits to an exponential function. For in- 
stance, for E = 0.3Vr (triangles), the fit gives Pe{J?) — 
(0.35/ctr) exp(-Q.12 3f/cr R ) for J*f > 30<T fl . (b) Average lo- 
calization length as a function of E c — E. The line is a fit to 
Eq. |6| where a, E c and v c are the fitting parameters. The nu- 
merical results (dots) have been obtained after averaging over 
fOOOO realizations of a blue-detuned isotropic speckle poten- 
tial in a box of length L — 800ct.r with grid step <tr/8. 



box center [if V(r = 0) < E] and Ae (r) = elsewhere. 
The calculation of the localization length, Eq. (J5|, at 
fixed energy is repeated for several realizations of the 
disordered potential. The corresponding probability dis- 
tribution Pg(_£?) is shown in Fig. 4la), for various en- 
ergies. It decreases exponentially for sufficiently large 
values of _£f. The average localization length, (Jjf(E)) = 
f d^f ££Pe is shown in Fig. El b) . It diverges when 
approaching the percolation threshold E c and, close to 
this value, it can be well fitted to the function 

(jSf(S))=a(fc^) "\ (6) 

where a = 1.63±0.3tr^, the critical exponent v c = 1.53± 
0.2 and the percolation threshold E c — 0.54 ± 0.02 are 
fitting parameters. The value of E c obtained from the 
fit is in agreement, within the numerical error, with the 
value E c = 0.52V R found with other methods [see Fig. [3]. 

Above the percolation threshold, for E > E c , not all 
the allowed energy regions extend to infinity and some 
particles can still be classically localized, depending on 
their initial position. This behavior is relevant to our 
system, where the dynamics starts inside the disordered 
potential. We have thus studied, as a function of the 
particle energy, the probability Pr to find, for different 
realizations of the disordered potential, an allowed region 
connecting two points at a distance R = \/ x 1 + y 2 . The 
results are shown in Figs. p[a) and (b) for the blue- and 
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Figure 5: (color online) Probability to find a classically- 
allowed region stretching to a distance larger than R. Each 
point has been obtained after averaging over 500 realizations 
of the disordered potential. Different lines refer to different 
values of R (indicated in the figure). Panels (a) and (b) cor- 
respond to blue- and red-detuned speckle potentials, respec- 
tively. The vertical dashed lines are the percolation threshold 
E c = 0.52V„. 

red-detuned speckle potentials, respectively. For large 
enough distance i?, the probability shows a sharp transi- 
tion, from zero, for E < E c , to a finite value, for E > E c . 
Note that, for E > E c , the probability converges to a 
curve which does not depend on the length R. This is 
consistent with the existence of allowed regions of finite 
size above the percolation threshold. The number of such 
regions rapidly decreases when increasing the energy. In 
a blue-detuned speckle potential [see Fig.[5]ja)], there are 
always forbidden regions but, for energy E > V R , almost 
all percolation regions extend to infinity (Pr quickly ap- 
proaches the value 1). In a red-detuned speckle poten- 
tial [see Fig. ||b)] there is a single percolation region for 
E > 0, as there is no forbidden region in this case. 

III. TRANSPORT REGIMES IN A 2D 
DISORDERED POTENTIAL 

In this section, we study the dynamics of a classical 
particle of fixed total energy E in a 2D anisotropic dis- 
order. The central goal is to characterize the behavior of 
the mean square displacement as a function of the par- 
ticle energy. We mainly focus on blue-detuned speckle 
potentials but also discuss results for the red-detuned 
case. 



A. Classical equations of motion 

The dynamics of a classical particle in a disordered 
potential is most conveniently written by rescaling the 
dynamical variables in the following way: 



E = £/|V R |, (7a) 

f = r/cr R , (7b) 

P = p/Pr, (7c) 

t = t/r R , (7d) 
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where E is the particle energy, r = (x,y) its position, 
P = (PxiPy) its momentum, and t is the time. In- 
troducing the momentum p R = ymfT^J and the time 
r R = yj ma R /\ Vr\, where m is the atomic mass, the clas- 
sical equations of motion read: 



df 

dF = p ' 

dp dv{r) 
At ~ df 



(8a) 
(8b) 



where v(r) = V(r)l V R \ is the reduced disordered po- 
tential. Equations (|8) highlight the universality class of 
the classical dynamics: upon the above rescaling by the 
disorder parameters V R and <7 R , the classical dynamics 
only depends on the properties of the reduced potential, 
w(f), that is on the model of disorder, and in particu- 
lar on its anisotropy. It is then sufficient to study the 
particle dynamics for a single set of parameters V R and 
(Jr. This is an important difference compared to the 
quantum dynamics, which also depends on V R / E R , where 

e r = h 2 /mal |22l ESS ED]. 

The reduced equations (J8|, for the speckle potential 
described in Section [nj form the basis of the calculations 
reported in this paper. For each realization of the dis- 
ordered potential V(r), we numerically solve Eqs. ^ to 
calculate the classical phase-space trajectory (r(t),p(t)), 
depending on the initial conditions (ro,p ). We use an 
adaptive ordinary differential equation algorithm where 
the time step is adjusted such that the energy is con- 
served within about 0.5% over the entire evolution. Nu- 
merical simulations are run with the fixed initial posi- 
tion ro = 0. The initial momentum is chosen with am- 
plitude |p | = y/2m[E — V{tq)\ and, unless specified, 
along a random direction chosen homogeneously. Dis- 
ordered speckle potentials are randomly generated and 
disregarded if V(tq) > E. They are typically created 
in a box of linear length L — 400cr R with a grid step 
Ax — Ay — cr R /10 and repeated periodically over the 
infinite plane [88] . Quantities of interest (see below) are 
then obtained by averaging over different realizations of 
the disordered potential. 



B. Time scale of ballistic dynamics 

The primary effect of the disordered potential is to 
modify the ballistic dynamics of particles. A particle of 
energy E initially moving along the r G {x, y} direc- 
tion looses the memory of its initial momentum within a 
characteristic time scale r* (E) [5T] . In order to identify 
r* (E), we have studied the momentum covariance tensor 



Cov r y(t) = (p r (t)p r >{0) ). 



(9) 



The typical behavior of the diagonal elements of 
the normalized momentum covariance, Cov r (t) = 
Cov r)I .(t)/(^(0) ), is shown in Fig.^a) as a function of 
time and along the r = x,y directions. We have checked 



that the off-diagonal elements of Eq. ^ vanish, in our 
case. The quantity Cov r (t) is characterized by a rapid 
decay and we define t*(E) such that Cov r [r*(-B)] = 0.1. 
For an anisotropic speckle disorder, we find t*(E) ^ 
t*(E) and t*(E) increases with oy (r g {x,y}). We 
then introduce t*(E) = min[r*(E),T*(E)]. For times 
t < t*(E), the dynamics is ballistic and isotropic. For 
times t ~ t*(E), the dynamics is modified by the disor- 
der, it enters various transport regimes and may become 
anisotropic (see Sec. HID). 



Note that, in our simulations, Cov r (t) does not ex- 
actly converge to zero at large times as exemplified in 
Fig. [6]ja). We indeed find fluctuations with both positive 
and negative values, which survive massive configuration 
averaging. We have however checked, from our numerical 
results, that the integral l r (t) = f Q du Gov r ^ r (u) /m 2 con- 
verges in the long time limit to a finite value (up to fluctu- 
ations), as shown in Fig. |6^b). This holds for particles of 
energy both above and below the percolation threshold. 
In addition, for particles in the normal diffusion regime, 
the limit of this integral equals the diffusion coefficient 
obtained with the fit of the mean-square displacement as 
discussed in Sec. Ill D[ ) . This is consistent with the cele- 



brated relation, D r y = J Q °° dt Cav r y (t) / m 2 , valid for 



normal diffusion |84) , where D r r , 
of the diffusion tensor. 



is the (r, r') component 



C. Energy redistribution 

Another main effect of disorder is to redistribute, after 
a time of the order of t*(E), the different forms (kinetic 
and potential) of the particle energy, under the constraint 
that the total energy is conserved. An example is shown 
in Fig. [6^c) where we plot the average kinetic energy 
along the x direction, (p x (t) 2 ) /2m, as a function of time. 
Different lines correspond to particles with the same en- 
ergy E = V R and different initial momentum directions 
p . After times of the order of t*(E), the lines converge 
to the same value, (E^), which will be calculated below 
[see Eq. (fl2j)]. 



A first insight into the dynamics is obtained by us- 
ing the statistical micro-canonical ensemble, which as- 
sumes phase-space ergodicity in the shell of energy E. 
The probability distribution of the potential energy E p 
for a particle of energy E in a single realization of the dis- 
ordered potential is then n(E p ) cx J dr dp S[E — p 2 /2m — 
V(r)] <5[V^(r) — E p ]. In a continuous, homogeneous disor- 
der, averaging over realizations of the disordered poten- 
tial equals spatial averaging, so that J dr (5[V(r) — V] oc 

p(v) m\- 

For a blue-detuned speckle potential, using Eq. ([2| , we 
find, 



n(E p ) 



e -E p /V R 



V R (1- e- E / v *) 



Q(E P )Q(E - E p ), V R >0, (10) 



normalized such that J dE p n(E p ) = 1. The distribu- 
tion n(E p ) is shown in Fig. ffid) for E = V R > 0. The 
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Figure 6: (color online) Primary effects of the disordered po- 
tential on the dynamics of a classical particle: modification 
of the ballistic dynamics and redistribution of energy. Here 
we consider particles of energy E — V R . Panel (a) shows the 
normalized momentum covariance Cov r (i) for particles with 
momentum p r (0) = \j2m[E — V(ro)] along the r — x (blue 
line) and r = y (red line) direction, as a function of time, (b) 
Integral of the momentum covariance, I r (i), (in units o~ r /t r ), 
as a function of time and along the r = x (blue line) and r = y 
(red line) direction, (c) Average kinetic energy, (p x {t) 2 ) /2m. 
The different solid lines refer to particles with fixed energy 
and di ffere nt initial momentum directions. The dotted line 
is Eq. (121. Analogous results are obtained for (p y (t) 2 } /2m. 



(d) Potential energy distribution obtained at time t = 30r R 
The dots are numerical results and the solid line is Eq. ( 10 1 



The numerical results have been obtained for a blue-detuned 
speckle potential (Vr > 0) by averaging over 10000 realiza- 
tions in panels (a) and (d) and 2000 realizations in panels (b) 
and (c). 



good agreement between the numerical data (dots) and 
the analytic prediction (Eq. (10); solid line) legitimates 



the micro-canonical statistical model above, and, con- 
sequently, the phase-space ergodicity hypothesis. Using 
Eq. (10), we find that the average potential energy of a 



particle of energy E in the disordered potential is 

E/V R 



(E p ) =V K [1 



e E/V R _ l 



Vr >0. 



(11) 



Then, using the relation imposed by energy conservation, 
E = (Ek) + (E p ), we find that the average kinetic energy 
is given by 



(£ k ) = V R 



E/V R 



1 - e-E/Vn 



1 , v R > 0. 



(12) 



In particular, in the limit E <C V R , we have (E p ) ~ E/2 
and (E'k) — E/2. This is easily interpreted: The particles 
are trapped in minima of the disordered potential, which 
can be assimilated to local (2D) harmonic oscillators, and 
the total energy is equally partitioned over the four de- 
grees of freedom. In the opposite limit, E ^> V R , the 
energy of the particles exceeds the disordered potential, 



the velocity field is nearly uniform and the average po- 
tential energy equals the spatial average of the disordered 
potential. We have (E p ) ~ V R and (E k ) ~ E - V R . Fi- 
nally, Eq. ( 12 ) is in excellent agreement with the asymp- 



totic value of the kinetic energy obtained in numerical 
calculations [see Fig. [fnc)]. For a red-detuned speckle 
potential, we find a similar behavior. The counter-parts 
of Eqs. (flOl), (fill and Jl2| have been calculated [33] and 



checked numerically. 



D. Energy-dependent transport regimes 

In the following we study the mean-square displace- 
ment, (r(t) 2 ), along the r £ {x, y} direction. Note 
that, after disorder averaging, the two parity symmetries 
x O —x and y O — y of the 2D (anisotropic) speckle 
potential impose (x(t)) = (y(t)) = (x(t)y(t)) = 0, for 
( > t' (E) when fixing the initial momentum direction. 
The above mean values hold at any time when averaging 
over an isotropic initial momentum distribution. 

Figures Jljd)-(f) show the typical behavior of (r(t) 2 ') 
as a function of time, for three different energies £ in a 
blue-detuned speckle potential (V R > 0). Quite generally, 
the plots in log-log scale show straight lines for long time 
periods. In each of them, the dynamics is well described 
by the general formula 



{r{tf) 



= 2D r x 



(13) 



where the parameters D r and a r depend on the rescaled 
energy E/V R . The different transport regimes are deter- 
mined by the value of the exponent a r : (i) classical local- 
ization (a r — 0), where (r(t) 2 } is independent of time; 
(ii) sub-diffusion (0 < a r < 1); (iii) normal diffusion 
(a r = 1); (iv) super-diffusion (1 < a r < 2); (v) ballistic 
expansion (a r — 2), which holds for free particles. 

As discussed above, for times t < t*(E), the dynamics 
is isotropic and ballistic. We have 



(14) 



where v = y/2(E^) a /\V R \ and (-Ek)o is the kinetic energy 
averaged over an isotropic momentum distribution. For 
the 2D geometry considered here we have (-Ek)o = (-^k)j 
which for a blue-detuned speckle potential is given by 
Eq. |l2]) [94 . Equation Q is the dotted black line 
in Figs, [ljd)-(f) and, as expected, it shows a very good 
agreement with the numerics [94j for short times. 

For longer times (t 3> t* (E)), the dynamics is strongly 
affected by the disorder and we have a r < 2. In the 
following, we discuss the transport regimes accessible to 
the particles, depending on their energy. 

Classical localization regime. Let us first discuss 
the dynamical behavior of a particle of energy below the 
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Figure 7: (color online) Integrated, average spatial distribu- 
tion (dots) along the x (upper row) and y (lower row) di- 
rections in semi-log scale. Different sets of data in each fig- 
ure refer to different evolution times, vertically shifted for 
clarity and explicitly indicated. Here we consider a blue- 
detuned speckle potential. In panels (a) and (b) the en- 
ergy is E — 0.3Vr. The solid lines are fit to the function 
f r (r,t) w e _iv ^' r ' which reproduces quite well the numerical 
results in the tails of the distributions. In particular, we have 
b x (t/r R = 10 4 ) = 0.11/crjj and b y (t/T„ = 10 4 ) = 0.22/a R , the 
ratio between the two values being equal to the anisotropy 
factor A = 2. In panels (c) and (d) the energy is E — V R . 
The dotted lines are guides to the eye. The black lines are 
fits with exponential functions, f(r,t) « e _i> *" ( - t ^ r ' , for rel- 
atively short times (t/r R = 150), while the dashed lines are 
f(r, t) « e~ r / 4D r( E )t^ reproducing the data for long times. In 
panels (e) and (f) the energy is E — 3Vr and the dashed lines 
are f(r,t) « e -r 2 /4Dr(E)t_ The diffusion coemc i e nts D r (E) 



are given by Eqs. ( 15 1 and ( 16 1 



percolation threshold (E < E c ). Figure |T|(d) shows that, 
in average, {r(t) 2 ) is well described by Eq. (13), with 
different parameters in different time intervals. At in- 
termediate times, the particle chaotically explores the 
finite-size allowed region it belongs to, bouncing between 
the borders. Its probability distribution progressively 
fills this region and we then find a subdiffusive behav- 
ior (a r < 1). At asymptotically large times, when the 
probability distribution has filled all the allowed region, 
we find a r = 0, corresponding to cla ssica l localization. In 
this case, the coefficients D r in Eq. (13) are given by the 



average mean-square size of the allowed regions along the 
x and y directions [these arc plotted as horizontal solid 
lines in Fig. [ljd)]. In an anisotropic disorder, we find 
that (x 2 )/(y/ — A 2 , which is consistent with the fact 
that the allowed regions are anisotropic, with the same 
anisotropy factor as the potential. 

In order to get further insight on the relation be- 
tween this behavior and the topography of the allowed 
regions, we have studied the "probability of diffusion", 
/(r, £|ro, to), which is the probability distribution to find 
a particle at position r at time t assuming that it is 
at position r at time t . Due to time translation in- 
variance and space transition invariance (after averag- 



Figure 8: (color online) Diffusion coefficients D r (E) (dots; 
in units <7r/tr) along the r = x,y directions. The results 
have been obtained numerically from fits of (r(t) 2 ), accord- 
ing to Eq. (13 1, for large times. The main figure shows the 
results for a blu e detuned speckle. Lines are t he p ower-law 
fit with Eq. (1151) and coefficients given by Eq. (161. The in- 



set shows the diffusion coefficients D r (E) (dots) for a red- 
detuned speckle disorder, as a function of _E/|Vk| + 2. Lines 
are a power-law fit with Eq. ( |17[ ) and coefficients given by 
Eq. (181. Numerical results have been obtained after averag- 



ing over 2000 realizations of the disordered potential. 



ing over disorder configurations), it is sufficient to con- 
sider the case i"o = and to = 0. Numerical results 
for the reduced probabilities f x (x,t) — J dyf(r,t\r ,to) 
and f v (y,t) = J dxf(r,t\r Q ,t ) are shown in Fig. [Wa) 
and (b) for particles of energy E — 0.3Vr in a bTue- 
detuned speckle potential. The distributions along x and 
y present tails that can be well fitted by exponential func- 
tions. This behavior can be quantitatively related to the 
probability distribution of the sizes Jz? of the allowed re- 
gions for a given energy E, which is exponentially small, 
Pe{-&) oc exp(— a E Ji?/<j R ) (see Sec. IIB). Consider par- 
ticles trapped in lakes (allowed region) of linear length J>f 
along x. They chaotically explore the lakes and the lin- 
ear density along x (averaged over the disorder) is almost 
uniform and scales as 1/Jzf. Hence, the "probability of 
diffusion" scales as f x (x,t — > oo) ~ J°° dS£ Pg(^f)^. 
We then find that, for x ^> a K /a,E, fx(x,t — > oo) de- 
cays exponentially with the same characteristic length 
of Pe{^), namely ln[f x (x,t — > oo)] ~ —a,Ex/<j R . For 
E = 0.3Vft, we indeed find that, in the long time and long 
distance limits, the decay rate of f x (x,t) [b x ~ 0.11/<7 R ; 
see Fig. [7ja)] is approximately the same as that of Pe 
[<3,e — 0.12/<t r ; see the red triangles in Fig. |4f a)] . The 
decay of f y (y, t) [see Fig. ^b)] is simply related to that 
of f x {x,t) by taking into account the anisotropy fac- 



tor A = (T x /a y of the disorder. Therefore, the function 
f v (y, t) has exponential tails which decay A times faster 
than f x (x,i). 



Normal diffusion regime. For sufficiently long 
times and for energies above the percolation threshold, 
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the dynamics is diffusive (a r = 1). This is shown in 
Fig. Ill e) and (f) where the solid lines are fits of (r(t) 2 ) = 
2D r \~E/V R ) x t to the numerical data for t t* , with 
D r as fitting parameter (along the r — x,y directions). 
We have studied the behavior of the diffusion coefficients 
D r (E/V K ) as a function of the particle energy, for blue- 
and red-detuned speckle potentials. 

For a blue-detuned speckle potential, we find a clear 
power-law scaling of the diffusion coefficients versus the 
energy, 



D r (E/V K ) = D° 



K. > 0, 



(15) 



as evidenced by Fig. [8[ w hich shows D r versus E in log- 
log scale. Fits of Eq. ( 15 1 to the numerical data, with D° 



and 7 r as fitting parameters provide 

D° x ~ (0.81 ± 0.03) <7 b /t R) V R > 0, 
D° y ~ (0.26 ± 0.01) a R /r R , 



1x — 7y 



3.15 ± 0.04. 



(16a) 
(16b) 
(16c) 



The diffusion coefficients D r are characterized, within 
numerical error, by the same exponent 7 along the x and 
y directions. In particular, the anisotropy of the diffusion 
does not depend much on the energy |95j . Interestingly, 
we find that, in the diffusion regime, the anisotropy factor 
of the diffusion, D x /Dy ~ 3.1, significantly differs from 
the squared anisotropy factor of the disordered potential 
(A 2 = 4). It shows that, differently from the classically 
localized regime, the anisotropy of the diffusion differs 
from the anisotropy of the disorder correlation function. 
This is further investigated in Fig. [9j which shows the 
ratio (D x /D y )/X 2 as a function of A = a x /a y for E = 
3V R . For A > 1, we find that D x (E)/D y (E) < A 2 . 
For a red-detuned speckle potential, we have 



D r (E/V R ) = D° r 



E-2V a 



Vr < 0. 



(17) 



Fits of Eq. ( 17 ) to the numerics with D® and j r as fitting 
parameters are shown in the inset of Fig. [8| We find 



D 

D° 

y 



~ (1.19 ± 0.07) a R /r R , V R < 0, 
- (0.37 ± 0.02) a R /r R , 
7^-7^- 3.45 ± 0.09. 



(18a) 
(18b) 
(18c) 



Compared to the blue-detuned case, we obtain signif- 
icantly larger coefficients D® and larger exponents j r . 
However, the anisotropy factor is approximately un- 
changed, D°/D° ~ 3.2. 

The power law scaling of the diffusion coefficients ver- 
sus the particle energy in large energy windows is a rather 
general feature of particle diffusion problems [IHl [Ml EI] • 
Here, we find that, remarkably, the exponents j x y are 
about the same along the two directions. For the energy 
window we are considering here, which is relevant to the 
experiment of Ref. [77], the disorder is strong. This is 




Figure 9: Anisotropy of the diffusion coefficient. Ratio 
[D x /D y ]/X 2 as a function of A = a x /a y . Diffusion coeffici ents 
are obtained from the fit of (r(t) 2 ), according to Eq. (131 
and for asymptotically large times. Dots are obtained tor 
E — 3V R and by averaging over 2000 realizations of a blue- 
detuned speckle potential. The dashed lines is a guide to the 
eye. 



suggested by the fact that the diffusion coefficients at 
a given energy are significantly different for blue- and 
red-detuned speckle potentials. For higher energies, the 
theory of quantum diffusion shows that the diffusion co- 
efficients, calculated in the Born approximation, only de- 
pend on the autocorrelation function of the disordered 
potential [69], which is the same for the blue- and red- 
detuned speckle potentials used in the present study. For 
isotropic Gaussian correlation functions, one then finds 
7 = 2.5 EH [75]. Note however that the latter result 
holds for E 3> V 2 /E B , which can be equivalently writ- 
ten as (E/V R ) 2 ^> (ka x>y ) 2 . The parameters used in the 
present work do not fulfill this condition since the valid- 
ity of Newton's equation of motion requires ka x ,y ^> 1 
and our calculations are limited to E < 10|V^|. 

In the normal diffusion regime, we expect that the 
probabilities of diffusion are given by Gaussian functions 
[6"?] . In Fig. 7[[e) and (f), we plot the (non-normalized) 
probability distributions along the x and y directions, 
taken at different times t = 400t r and t = 800r R , for 
a blue-detuned speckle potential and for particles ini- 
tially at the origin, ro = 0. The numerical findings 
are in very good agreement with the functions f r (r, t) cx 
e -r /4D r (E)t^ wnere coefficients D r {E) are given by 
Eqs. (15) and (16). We have checked that this holds for 



different energies in the normal diffusion regime. 

Finally, we have studied the correlation functions 
Ck(E,t) — (x k (t)y k (t)) / (x k (t)) (y k (t)) for various integer 
values of k. For particles in the diffusion regime and for 
t » t*(E), we find Ck(E,t) — > 1 for different values of 
k. This shows that, in the normal diffusion regime, the 
x and y variables are independent. Then, taking into ac- 
count the results of Figs, ^e) and (f), the probability of 
diffusion, i.e. the probability density to find a particle at 
coordinate (x, y) which is in Tq — (xo,yo) at time to, is 
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Figure 10: (color online) Scaling coeffi cien t a r (E) obtained 
from a fit of (r(i) 2 ) [according to Eq. (13)] at intermediate 
times, 100 < t/r R < 1000, (circles) and at larger times, t/r R > 
1000, (triangles). The fit are shown, for instance, in Fig.JTJd)- 
(f). Solid and dotted lines are guides to the eye. Different 
colors refer to r = x (blue lines and filled symbols) and r = y 
(red lines and empty symbols) direction. Panels (a) and (b) 
refer to blue- and red-detuned speckle potential, respectively. 
In each panel, the vertical dashed black line is the percolation 
threshold, E c = 0.52V R . 



given by 



f E {r,t\r ,t ) = 



(y-vo) 



g 4B x (i3)(t-t ) g 4D s (E)(t-t ) QU _ 



Sub-difFusion regime. Anomalous diffusion occurs 
in a variety of physical systems jHH H7J [35] and has been 
experimentally studied in Refs. [MJ |jJD - f5§] for instance. 
It is denned by the condition a r ^ 1 in Eq. ( 13 ) and 



two physically different regimes should be distinguished: 
sub-diffusion when < a r < 1 and super-diffusion when 
1 < a r < 2. In our case, super-diffusion only shows 
up for t ~ t* in the crossover from ballistic expansion 
(a,. = 2 for t < t*) to normal or sub-diffusion (for 
t ^ T *)- Sub-diffusion is more relevant to our study. 
For a blue-detuned speckle potential, a transient sub- 
diffusive behavior is found for E < 2V R and long (al- 
though not asymptotic) time scales with i » t*. For 
instance, the fits of Eq. (13) to (r(t) 2 ) in Figs.jljd) and 

(e), at intermediate times and with D r and a r as fit- 
ting parameters, (dashed lines) give a r < 1. The scaling 
coefficients a r , as obtained from such fits in two dif- 
ferent time windows, are plotted as a function of en- 
ergy in Fig. 10 for blue-detuned (left panel) and red- 
detuned (right panel) speckle potentials. The value of a r 
obtained from the fits in an intermediate time window 
(10 2 t r < t < 10 3 r R ) smoothly crosses over from a r = to 
ot r = 1 when crossing the percolation threshold E c . For 



very long times (i > 10 3 t r ), the crossover is very sharp, 
which is compatible with an abrupt transition between 
two asymptotically-relevant dynamical regimes: Below 
the percolation threshold (E < E c ), the mean-square 
displacement (r(i) 2 ) evolves from subdiffusion to classi- 
cal localization [see Fig. [l|e)]. Conversely, slightly above 
the percolation threshold (E > E c ), (r(t) 2 ) evolves from 
subdiffusion to normal diffusion [see Fig. nlf)]. Notably, 
this case shows that the sub-diffusion regime cannot al- 
ways be considered a precursor of localization. 

While the subdiffusive behavior is not asymptotic, it 
is however experimentally relevant. In fact, subdiffu- 
sion is observed for E ~ V R , at times of the order of 
~ 200t r , which corresponds to ~ 100ms for the param- 
eters of Ref. [77], that is the typical time scale of the 
experiment. In the sub-diffusion regime, the probabil- 
ity density /(r, t|ro, to) is characterized by exponentially 
decaying tails [see Figs.jTja), (b), (c) and (d)]. We inter- 
pret this result by analogy with the behavior described 
above in the classical localization regime. Slightly above 
the percolation threshold, there exist allowed regions of 
finite size that are disconnected from each other or con- 
nected by very thin bottlenecks to a percolating cluster. 
Then, as illustrated in Fig. [Tl|a), a particle stays long 
times in the same lake where it behaves similarly as be- 
low the percolation threshold. In particular, its dynam- 
ics is subdiffusive. As discussed above, this persists for 
asymptotically large times when classical localization oc- 
curs [for E < E c ; Figs. [W a) and (b)]. Above the percola- 
tion threshold (E > E^) however, the particle can find a 
path connecting the lake to another one after sufficiently 
long exploration of the lake. It then spends a long time 
in the new lake before switching to another one, and so 
on [100) . For longer time scale and larger space scale, the 
dynamics results from random jumps of the particle from 
a trapping region to another one, hence developing nor- 
mal diffusion, similarly as in usual Brownian motion. In 
particular, the probability of diffusion becomes Gaussian 
[see Figs.[7|c) and (d) for sufficient long time]. Here, the 
existence of long waiting times would explain the tran- 
sient sub-diffusive regime obtained at intermediate times. 
Note that for higher energies, where no transient sub- 
diffusion is found, the trajectory of a particle does not 
show long waiting times and reminds that of Brownian 
motion even on a short times scale [see Fig. [TT|(b)] . 



IV. EXPANSION OF AN ATOMIC CLOUD 



Section |III| presents a detailed analysis of the trans- 
port regimes for particles with a fixed energy. In typical 
experiments with cold atoms, transport is rather investi- 
gated from the behavior of a cloud of atoms with a broad 
energy distribution. In this section, we study the expan- 
sion of an atomic gas by considering the atoms as non- 
interacting classical particles moving in a 2D disordered 
potential. We take into account the energy distribution 
of the gas and discuss our findings in the light of the 
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Figure 11: Trajectory of a particle of energy E in an 
anisotropic, blue-detuned speckle potential, obtained by solv- 
ing Eqs. (pp. In the panel (a) the energy is E = V R and the 
trajectory is characterized by relatively small regions where 
the particle spends long times. This behavior give rise to 
transient sub-diffusion. The panel (b) shows the trajectory 
for a particle of energy E — 3Vr. The dynamics is diffusive 
and qualitatively similar to a Brownian motion |67| . 



results of Sec. HIT 



A. Energy distribution of an initially trapped gas 

The aim of this section is to calculate the energy dis- 
tribution of a thermal gas of non-interacting particles 
initially trapped. We will consider two experimentally 
relevant cases in which the disorder is initially combined 
with the trap or switched on after releasing the trap. The 
energy distribution is highly sensitive to the way the dis- 
order is turned on. 

Quite generally, in transport experiments with cold 
atoms, the atomic cloud is initially in equilibrium in a 
trap described by the potential Vb(r). The initial single- 
particle classical Hanhltonian is Ho(r, p) = |p| 2 /2m + 
Vq (r) and the phase-space density distribution for a ther- 
mal gas at temperature T reads 



Po(r,p) = 



~Zo 



(20) 



where Z = J dr'dp' e -PH (r',p') ^ p = i/^T and k B 
is the Boltzmann constant. At time t = 0, the con- 
finement is switched off and the non-interacting gas ex- 
pands in the disordered potential V(r). The dynam- 
ics is then governed by the single-particle Hamiltonian 
H(r,p) — \p\ 2 /2m + V(r) and the joint distribution of 
energy and position at time t = + is 



n(v,E) = J dp Po (r,p)6(E-H(r,p)) (21) 

*-e{E-V(r)). (22) 



e -P(,E-V(r)+V (r)) 



rdr'e-W') 

In the following, we calculate the disorder-averaged en- 
ergy distribution (n(r, E)) 1 distinguishing the two exper- 
imentally relevant situations: (i) The gas is prepared in 
a bare trap, so that Voir) = V^ rap (r), where ]/ trap (r) is the 



potential of the trap, which is independent of the disor- 
dered potential. The disorder is suddenly switched on 
after the gas has thermalized in the trap, (ii) The gas is 
prepared in the trap and in the presence of disorder, so 
that Vb(r) = V tlap (r) + V(r), which now depends on the 
disordered potential. 

Case (i). If Vo(r) does not depend on the disordered 
potential, averaging the distribution Eq. (22 1 gives 



i(v,E)) = 



(3e 



-/3V (r) 



J dr / e -pv a {V) J 



dVP(V)r 



-P(E-V) 



where P(V) is the one-point probability distribution of 
the disordered potential. It is interesting to notice that, 
in a homogeneous disorder and independently of the spe- 
cific form of the initial trapping potential, the energy and 
position variables decouple: 



(n(r,E)) = n tiap (r)(n(£)), 



where 



n tiap (r) = / dpp (r,p) 



-PV (r) 



/di-'e-W') 
is the initial spatial distribution and 



i(E)) =0 f dVP(V), 

J — OO 



-P(E-V) 



(23) 
(24) 

(25) 



is the average energy distribution in the presence of dis- 
order. For a speckle potential, Eq. ^ holds and the 
integral in Eq. (25) can be easily done. In particular, for 
a blue-detuned speckle potential, we find (104] 



(n(E)) = 



P 



1-/3V R 



e -PE _ e -E/V R 



Q(E). 



(26) 



A plot of Eq. ( 26 1 as a function of the energy is shown 



in Fig. 12 [a,) for different values of f3V R . 



Case (ii). If on the contrary, Vb(r) depends on the 
disordered potential, the calculation of Eq. (22 1 is diffi- 



cult since both the numerator and the denominator de- 
pend on the specific disorder configuration. A simpli- 
fied expression is obtained if the thermal gas extends, in 
the initial trapping potential Vb(r), over a region much 
larger than the correlation length of the disorder. This 
condition can be fulfilled for a sufficiently shallow con- 
fining trap. In this case, the self-averaging properties 
of the disordered potential allow us to substitute the 
integral J dr' e _/3y °^ r ' by the disorder-averaged value 
(/ dr' e~W r ')) = / dr'e~^ r ') / dVP(V) e -? v . Wc 
thus find that, also in this case, the energy and posi- 
tion variables decouple. Then, Eq. (231 still holds, with 
the same density distribution, given by Eq. ( 24 1 , but a 
different energy distribution, 



(n(E)) 



(3e^ E J^dVPjV) 



(27) 
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Figure 12: Disorder-averaged energy distributions of the 
atomic cloud. Panel (a) is Eq. ( |26[ ), obtained when the speckle 
potential is suddenly switch on after releasing the trap where 
the atomic cloud is initially in equilibrium [Case (i)]. Panel 
(b) is Eq. ( 28 1 and corresponds to a thermal cloud in equilib- 
rium in a disordered trap [Case (ii)]. The different curves refer 
to different values of V R : Vr = 5ksT (solid line), Vr = k&T 
(dashed line) and Vr = 0.2 fee T (dot-dashed line). The gray 
shaded region indicates sub-diffusive particles (E < 2 Vr), the 
dar kest region classically localized particles (E < 0.52 Vr) [see 
Sec.~ 



III 



Inset: density of particles with energy smaller than 
E (solid lines). The dashed lines are a prediction in the large 
temperature limit k^T ^> Vr, see footnote |105| . 



In particular, for a blue-detuned speckle potential, we 

"l _ e -E/V R 



have [104] 



Q(E). (28) 



In the case of an isotropic harmonic trap K ra , p (r) 



|r| 2 /2, we find (r) 



T /27T(Ty, where 



<?t = V kftT j 'mui 2 . The self-averaging condition invoked 
above is thus cry >(Tjt- A plot of Eq. (28 1 as a function of 
energy is shown in Figjl2jb) for different values of (3V R . 
Equations (26) and (28 1 tend, in the high-temperature 
limit, k B T > E, Vr, to[n(E)) = 0(1 - e~ E ^). At low 
temperature, k B T <C V R , the two distributions differ sub- 



stantially, Eq. ( 28 1 being more strongly peaked at small 



energies. 

The insets of Fig. [12] show the fraction of particles 
with energy lower than E, n < (E, T) = J Q de (n(e)} (105] . 
as a function of temperature and for various energies. In 
particular it shows the fraction of atoms in the classically 
localized regime, n< (E — 0.52V R , T), and below the tran- 
sient sub-diffusive regime threshold, n < (E = 2V R ,T), 
identified in Sec. |III| The dotted line is the asymptotic 
behavior for large temperatures, n < (E,T) = (e~ E / Vli + 
E/V R — 1)/3Vr, which is the same for both cases (i) and 
(ii). When the gas thermalizes in the trap and in the 
presence of disorder [case (ii) described by Eq. (28)], the 



case when the disordered is suddenly ramped after releas- 
ing the confining trap [case (i) described by Eq. (|26[)]. In 
case (ii), all particles can be classically localized at low 
temperature [i.e. n < (E = 0.52Vr,T) ->■ 1 when T 0]. 
Conversely, in case (i), the speckle potential is switched 
on after releasing the initial confining trap, which trans- 
fers an average energy ~ V R to all particles and the frac- 
tion of particles at low energy is depleted. 



B. Expansion of the atomic gas 

We now study the dynamics of an atomic cloud ini- 
tially trapped in a pure harmonic potential [i.e. without 
disorder, case (i)]. At t — the trap is released and the 
atoms expand in a blue-detuned speckle potential. As 
shown above, the energy and the initial position of the 
atoms are uncorrelated variables, when averaging over 
disorder configurations. In our simulations, the energy 
is randomly chosen according to Eq. (26), while the ini- 



tial position is randomly chosen according to a Gaussian 
distribution of width gt = 10er R , which is relevant to 
typical experiments [77] . We expect this approach to be 
reproducible in a single run of the experiment since, for 
ctt ~> <7r, the speckle disorder "self-averages" over the 
atomic cloud. 

Density distribution. An important feature of the 
experiments with cold atoms is the possibility to ob- 
tain an image of the density cloud at different expan- 
sion times. For instance, it has been crucial in the recent 
observation of Anderson localization with Bose-Einstein 
condensates in the quasi- ID geometry [6"T1 RS] . 

The time-dependent density profile (averaged over dif- 
ferent realizations of the disordered potential) is calcu- 
lated by integrating the probability of diffusion over all 
energies and all initial positions 



n(r,t) = J dr / d-En tr 



,(r o )(n(E))f E (r,t\r o ,0), 



where n trap (ro) is the spatial density distribution in the 
initial trap. As noticed in the previous paragraph, for 
a thermal gas, the initial joint position-energy distribu- 
tion decouples into the product of separate distributions. 
This feature is particular to the case of an initial ther- 
mal distribution. For instance, it does not hold for ini- 
tially expanding Bose-Einstein condensates as considered 
in Refs. [3TJ [35J [7S] . The quantity /^(r, t\r , 0) gives the 
probability distribution to find a particle of energy E at 
position r and at time t > 0, assuming it was in ro at time 
t = 0. Its behavior has been discussed in Sec. |III| for par- 
ticles of fixed energy in various transport regimes. The 
integrated density profiles at time t = 400t r are shown as 



fraction of particles at low energy is larger than in the 



thick solid lines in Fig. 13 It is striking that these profiles 
are in general, not Gaussian. They present a pronounced 
central peak, which is partially due to the contributions 
of the atoms in the classical localization and sub-diffusive 
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Figure 13: (color online) Integrated spatial density distribu- 
tion of an expanding thermal gas in a 2D disordered potential. 
The thick solid line is the integrated spatial density distribu- 
tion (a) n x (x,t) = j Ayn(r,t) and (b) n y (y,t) = J dxn(r,t), 
where n(r, t) is given by Eq. p9| ). Here, V R = k B T, t = 400r R , 
the initial distribution is a Gaussian of w idth ctt = 10or and 
the energy distribution is given by Eq. ( 26 1 . In the numer- 
ical simulation, we used 10000 particles. We also plot the 
contributions of particles with energy E > 2V R (dashed line) 
and with energy E < 0.52Vr (thin solid line). Note that for 
the time considered here, the transient sub-diffusive regime is 
relevant and the fraction of atoms in this regime cannot be 
neglected. 



400 



D 




Figure 14: (color online) Mean square size, [Ar(i)] 2 , as a 
function of time, in the r = x (solid blue line) an d r = y 
(solid red line) directions. The dotted line is Eq. (301 and has 



been obtained by assuming a ballistic expansion for all atoms. 
The dashed lines are Eq. (310 , corresponding to all particles 
behaving diffusively. The horizontal line is the mean square 
size of the initial cloud, given by a\. Here V R = k B T and 

<7T = W<TR. 



regimes. It is worth stressing that, even taking into ac- 
count only the fraction of atoms in the diffusive regime, 
for which the probability of diffusion, /js(r,t|ro,0), is a 
Gaussian function [see Eq. (19)], the overall distribution 
is, in general, not Gaussian. This is due to the integration 
over the energy variable in Eq. ( 29 ) and to the strong en- 



ergy dependence of the diffusion coefficients [751 uM ■ In 
fact, since the particles in the diffusive regime determine 
the long tails of n(r,t), the spatial distribution can be 
used to extract the diffusion coefficients from experimen- 
tal data (see Ref. [77] and Sec. |v|. 

Mean square size. A simple quantity which can be 
extracted from the density profiles is the rms size of the 

atomic cloud, Ar(t) = yj J dr r 2 n(r, t), in the r = x,y 

direction and where n(r,t) is given by Eq. (26). The 



evolution of the rms sizes in the two directions are shown 
in Fig. [l4j for V B = k^T. At sufficiently short times, 
t <C t*(E), where E is the typical energy of the atomic 
cloud |107j . the rms size can be obtained by assuming 
that all the atoms behave ballistically |106] . In this case, 
we find 



[Ar(t)f = 4 



/3Vk ' 



(30) 



where <tt — \J k^T jniuS 1 is the rms size of the initial 
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atomic cloud. Equation (|30|) is shown by the dotted black 
line in Fig. 



For long times, t ^> t*(E), the rms size 
of the atomic cloud can be calculated by assuming that 
all the particles behave diffusively. This approximation 
is legitimate in the long-time limit if the fraction of clas- 
sically localized particles is negligible, which strongly de- 
pends on the ratio k B T/V R . In this case, the evaluation 



of Eq. (291 using Eqs. (26 1, (19) and (15) gives 



[Ar(t)] 2 ~ 2D r a 2 R 



[i-Q3Kr- +1 ]r( 7r + i) t 

(l-pV R )(/3V R pr Tr 



(31) 



where T(x) is the Gamma function. In Fig. 14 



wo 



show Eq. (31) calculated for the parameters D® and 7, 



given by Eqs. (16). We find a good agreement with 



the numerical findings for long times, despite the fact 
that the fraction of localized or sub-diffusive atoms is 
n< (E = 2V R ,T) ~ 70%. Hence, Eq. plj) could be used 



experimentally to extract the values of j r and D® charac- 
terizing the diffusive regime, and thus verify the results 
of Sec. [mDl 



V. EXPERIMENTAL RESULTS 

The study of classical transport of cold gases in a 2D 
disordered potential is experimentally feasible. A 2D 
disordered potential can be created experimentally by a 
speckle field, as discussed in Sec. |II A| The atoms can be 
confined to a plane by a tight transversal harmonic trap 
of frequency lo z . The dynamics in the transverse direc- 
tion can be neglected if k^T is much smaller than the 
vertical confinement energy Hlo z . 

In recent experiments |77j . we have studied the dy- 
namics of a cloud of 87 Rb atoms in a speckle disor- 
der. The geometry is not strictly 2D since kBT rj 
k B x 200(20) nKsa 6huj z , where w 2 /2tt = 680 Hz. There- 
fore, taking into account the harmonic confinement, the 
external potential is given by V(x, y, z) + mw 2 z 2 /2, with 
V(x,y, z) = Vi so {xsmQ — zcos9,y) where Vi so (u,v) is a 
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2D isotropic blue-detuned speckle potential with corre- 
lation length a = 0.8 /im and Gaussian autocorrelation 
function. In the experiment, the speckle disorder is tilted 
by 8 — 30° with respect to the horizontal plane. This 
makes the disorder anisotropic in the expansion plane, 
with correlation lengths a x — 1.6 /im and a y = 0.8 /xm, 
and anisotropy factor A = cr x /a y — 2. The speckle is 
assumed to be invariant along its propagation axis. This 
is a reasonable approximation since the longitudinal cor- 
relation length, ai on g — 9 fim is large compared to the 
vertical size of the cloud, (k^T '/muj 2 ) 1 ^ 2 ps 1 /im. 

Using a classical particle model is a reasonable ap- 
proximation since ka x ~ 10 and ka y 5, where 
k = y^mfcflT ' jh. In addition, the "correlation energy" 
E n = h 2 /ni(T R w k B x 2 nK is typically more than one 
order of magnitude smaller than the amplitude of the 
disorder V B = k B x 53 nK. Within these conditions, the 
localization length should be much larger than the system 
size (~ 1 mm in each direction) and one can expect that 
quantum interference effects can be neglected. Due to 
the tilting angle of the speckle pattern, which couples the 
in-plane dynamics to the dynamics in the confined direc- 
tion, the system is 3D rather than purely 2D. Therefore, 
the classical equations of motions Q must be adapted 
to take into account the vertical confinement ITT 



Remarkably, the results of this 3D model differ only 
slightly from the results presented above for the purely 
2D situation. The two asymptotic regimes (classical lo- 
calization and normal diffusion) are unchanged since the 
3D potential results from a cylindrical stretching of a 2D 
potential. We also find a transient sub-diffusion regime 
in approximately the same energy window and time scale 
as for the pure 2D case, 0.52U R < E < 2V B . For each 
regime, the dynamics is qualitatively similar for both the 
pure 2D and the trapped 3D cases. The diffusion coef- 
ficients have a power law scaling analogous to Eq. ( 15 ) 
|77| . The numerical simulations give diffusion parameters 
slightly different to the ones of Eq. ( 16 ): 



D° x = (1.20 ± 0.05) al/r K , 
D° y = (0.33 ± 0.02) a 2 Jr B , 
j x ~j y = 2.8 ± 0.06. 



(32a) 
(32b) 
(32c) 



For the parameters of the experiment, V R = k& X 53 nK, 
a E = <J X = 1-6 /im and t r = 0.71 ms, we obtain D° = 4.3 
/im 2 /ms and D® = 1.2 /im 2 /ms. 

In the experiment, the speckle disorder is switched on 
after releasing the initial confining trap in the (x,y) direc- 
tion, while the transverse confinement is kept during the 
whole expansion. We expect that a few collisions redis- 
tribute the energy between the atoms during the first 10 
ms of the expansion. The calculation of (n(E)) is mod- 
ified compared to Eq. (28) in order to include the con- 
tribution of the different vertical energy levels. Taking 
into account the quantization of the transverse harmonic 




Figure 15: Disorder-averaged energy distribution in the ex- 
periments of Ref. [77]. The solid line is Eq. ( |33[ ). Gray re- 
gions indicate different dynamical regimes: classical localiza- 
tion for E < 0.52Vr (darkest region), transient subdiffusion 
for E < 2V R (gray region), diffusion for 2V R < E < 5.5V R 
(white region) and transient ballistic expansion (see text for 
details) for E > 5.5Vr (light gray region). Here the ampli- 
tude of the disordered potential is Va = fe x 53 nK and the 
temperature T = 220 nK. 



oscillator, we obtain 



i(E)) 



oc e 



E/hu> z 

-E/k B r 

n=0 



^nhui z -E)/V B 



(33) 



which is plotted in Fig. [l5j For the considered exper- 
imental parameters, the fraction of classically localized 
atoms is about 0.4% of the total number of atoms while 
the fraction of subdiffusive atoms for the time scale of 
the experiment is about 6%. 

The goal of the experiment is to extract the energy- 
dependent diffusion coefficients in the 2D anisotropic 
speckle potential. They are experimentally measured by 
analyzing the complete 2D column density profile, ob- 
tained experimentally by fluorescence imaging along the 
vertical axis (see Fig. 16 ). In order to fit the experimental 
density, we use Eqs. ( |19| ) (with to — 0) and ( 29 ) convolved 
by a Gaussian function which takes into account the res- 
olution of the imaging system. Using the fact that, for 
an initial harmonic trap, n trap (ro) is a Gaussian function, 
we find 



n(x,y,t) = / dE(n(E)) 



2(2D H (E)t + <r^) 



2nty/(2D x (E)t + al){2D y {E)t + a 2 

(34) 



where (n(E)) is given by Eq. (33) and <tq = 15 /Ltm as 
determined experimentally. Equation (34) assumes that 
all the particles behave diffusively. This is a good ap- 
proximation since, for the experimental parameters, the 
fraction of classically localized and subdiffusive particles 
is small. Diffusion coefficients are extracted by using 
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the diffusion coefficients. Note that this is not the case 
after 50 ms as ballistic atoms are visible (see Fig. 2 in 
Ref. [77]). 

In addition to the investigation of the normal diffusion 
regime discussed above [77], we have recently performed 
experiments with a lower temperature of 30 nK. In this 
case, Eq. ( 33 ) leads to 8% of classically trapped atoms 
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and 55% of atoms in the subdiffusive regime. We observe 
little expansion of the cloud as a function of time. For 
times of several seconds we observe losses due to resid- 
ual heating. We are unable to quantitatively analyze the 
shape of the cloud because of the limited resolution of 
the imaging system and of the initial size of the cloud. 
The subdiffusive regime and the percolation transition 
presented in this manuscript thus remain to be experi- 
mentally studied in 2D cold atomic systems. 



Figure 16: (color online) Experimental two-dimensional col- 
umn density distribution after 200 ms (= 281r R ) expansion of 
the atomic cloud. Dots are cuts along x (left panel) and y 
(right panel) directions at cut positions explicitly indicated, 
artificially offset for clarity. Lines are the fitting function 
Eq. (34 1. Note that the spatial distribution is not Gaussian. 



the power law scaling predicted numerically D r (E) 
D° r x [E/Vr)^ (r = x,y). Figure 
2D density distribution is well reproc 



1G 



shows that the 



uced by the fitting 
function, Eq. (34 1, with D® and j r as fitting parameters. 



Experimentally, we find: 



D° x = 5.5 ±2.7 /im 2 /ms, 
D° = 1.6 ±0.8 /W/ms, 
j x ~ ly = 3.3. 



(35a) 
(35b) 
(35c) 



These values are, within experimental errors, in fair 
agreement with the classical theory [see Eqs. (32)]. In 



particular, the algebraic scaling inferred from the nu- 
merical simulations is consistent with the experimental 
findings. In addition, these results show the possibility 
to extract the single-particle diffusion coefficients from 
the dynamics of an atomic cloud. 

We finally notice that the high energy atoms are ballis- 
tic on the time scale of the experiment (see Fig. 15 ). The 



crossover between the ballistic and the diffusive behavior 
takes place when the Boltzmann time t-q (= AD(E) jv 2 , 
where v = \J E/m is the rms in-plane velocity) is of 
the order of the expansion time t = 200 ms. The same 
crossover between ballistic and diffusive behaviors is seen 
as a function of time in Fig. |ljc) for a given energy. 
For the experimental parameters, we estimate that below 
-Ebai-dif ~ 5.5 Vr, the atoms are diffusive (calculated in 
the x direction which gives the lower estimate of -Ebai-dif ) • 
Above this value, the atoms may be ballistic on the time 
scale of the experiment. In this case they however extend 
over a distance vt > 1 mm, larger than the spatial size 
used in Fig. [16] where we fit our diffusive distribution. 
Ballistic atoms thus do not affect the result of the fit for 



VI. CONCLUSIONS 

In this work, we have studied the transport properties 
of a cold atomic gas in a 2D anisotropic speckle potential, 
in the classical regime. We have shown that, upon proper 
rescaling, the atomic dynamics can be parametrized by 
the energy as the sole relevant quantity. Then, we have 
identified two asymptotically-relevant regimes, depend- 
ing on the particle energy: (i) For energy below the per- 
colation threshold, any particle is trapped in a finite- 
size region, so that (r 2 ) — > const (classical localization 
regime), along the r — x,y direction. The probability 
of diffusion (defined as the average density probability 
of finding an atom at a given distance from its origin 
point) develops exponentially decaying tails, which cor- 
respond to the probability distribution of the size of the 
classically-localized regions. Due to the pure topograph- 
ical nature of the latter, the probability of diffusion is 
anisotropic with the same anisotropy factor as the speckle 
potential, (ii) For energy above the percolation thresh- 
old, there appear infinite-size allowed regions in which 
the gas expands diffusively, (r 2 ) oc t (normal diffusion 
regime). The probability of diffusion is a Gaussian func- 
tion with anisotropy factor (along the x and y directions) 
significantly different from the anisotropic factor of the 
speckle potential. The diffusion coefficients along the two 
main directions of the speckle potential strongly depend 
on the particle energy. More precisely, we found a power- 
law scaling, D r ~ E 1 with 7 ~ 3.15 for a blue-detuned 
speckle potential and 7 ~ 3.45 for a red-detuned speckle 
potential. 

In addition, a transient sub-diffusion regime, where 
(r 2 ) oc t Ur with a r < 1, characterizes intermediate time 
scales, relevant to recent experiments |77| . This regime 
appears below and slightly above the percolation thresh- 
old. It is interpreted as due to the dynamics of ex- 
ploration of finite-size allowed regions which are discon- 
nected from each other or connected by narrow bottle- 
necks. For energy below the percolation threshold, the 
long-time dynamics crosses over to classical localization 
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(a r = 0). Conversely, above the percolation threshold, 
it crosses over to normal diffusion (a r = 1). This ex- 
ample shows that subdiffusion should not be regarded as 
precursor of classical localization. 

Using the above results as a groundwork, we have stud- 
ied the expansion dynamics of a cloud of cold atoms with 
broad energy distributions, as relevant to recent experi- 
ments. First, we have calculated the energy distribution 
of the cloud in the presence of disorder. Then, we have 
discussed the dynamics. The time-dependent profiles are 
in general anisotropic and non-Gaussian, with tails gov- 
erned by diffusive atoms. We have proposed a method to 
extract the energy-dependent diffusion coefficients from 
these tails, which involve many energy components. The 
experiments of Ref. |77j are realized in a quasi-2D ge- 
ometry, which does not significantly change the above- 
discussed dynamics. The energy-dependent diffusion co- 
efficients extracted from experimental data are in fair 
agreement with numerical calculations performed in the 
precise experimental geometry. 

In the future, it would be interesting to experimentally 
study the classical localization regime as discussed in the 
paper. In particular, working at lower temperature, it 
would be possible to create a gas of atoms with a sig- 
nificant fraction below the percolation threshold. Then, 
studying the fraction of classically localized atoms ver- 
sus temperature, for instance, would give direct access to 
the percolation transition and, in particular, to the crit- 



ical exponent. Another interesting direction would be to 
extend the study of transport in anisotropic 2D speckle 
potentials in the quantum regime, where traces of weak 
and strong localization effects can be searched for. 
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